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1. Coordinates . - The shell shall he considered as a three- 
dimensional continuous medium; for the coordinate surface, the mid- 
dle surface of the shell shall he assumed parallel to the hounding 
surfaces . Let a and p he the curvilinear orthogonal coordi- 
nates of this surface, coinciding with the lines of principal cur- 
vatures, and 7 the distance along the normal from the point 
(a,p) of the coordinate surface to any point (<x,p,7) of the 
shell (fig. 1). 


The square of the line element in spatial orthogonal curvi- 
linear coordinates is given hy the formula 



da 2 dp 2 dy 2 

h l 2 h 2 2 h 3 2 


( 1 . 1 ) 


where h-^ = h 1 (a,p,7), h 2 = h 2 (a,p,7), and h 3 = h 3 (a,p,7) are 
the so-called differential parameters of the first kind represent- 
ing for the chosen coordinates given functions of a,p,y. 


In the triorthogonal system of coordinates chosen as indicated; 


1 = A(l+k l7 ) 

h 2 = _L_ 

2 B(l+k 27 ) 

h 3 = ! > 


( 1 . 2 ) 


where A = A(a,p) and B = B(a,p) are the coefficients of the 
first quadratic form; and where kp = kp(a,p) and k2 = k 2 (a,p) 
are the principal curvatures of the coordinate surface on the lines 
corresponding to p = constant and a = constant, respectively. 


*"Osnovnye Differentsialnye Uravnenia Ohshche Teorii Uprugikh 
Oholochek." Prikladnaia Matematika I Mekhanika. Vol. 8, 1944, 
pp. 109-140. 
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The magnitudes 
by the relations 


x l> h 2> 


A, B, k 


l> 



and k 2 are connected 
- kjkgAB 


, Sh 2 1 1 SB 

1 Sa ~ A da 

. Sh l _ 1 Sa 

2 Sp “ B S3 


(1.3) 


J 


obtained from the equations of Lamb for the differential parameters 
hi, h2, and h.3 defined by equations (1.2) and from the equa- 
tions of Codazzi 




_S_ 

S3 


(k lA > . k 2 | 


J 


(1.4) 


From equations (1.2), (1.3), and (1.4) follow the equalities: 

"\ 


h c 


-1 

Sh 2 Sh x 
Sa S7 

-1 

Shi Sh- 


-1 


-1 


S3 


S7 


S%_ 

da 5 7 


sS 

S3S7 


(1.5) 
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2. Fundamental Equations of Three-Dimensional Problem of Theory 
of Elasticity . - The six components of the strain tensor of a dense 
medium are determined, in the system of coordinates assumed, by 
the equations 
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e__ 

COCO 


du 

hi «£ + h!h 2 
OX 



+ 



e pp = h c 


3 up 


3 ho -1 


dho " 1 


sr + 111112 ~sr~ + hs ~w~ “ r 


du^, 

®77 = 


11 6 ^>Fl (hA) 


" /<x P h 2 da 


> (2.1) 


l S / \ du> 

= JT s; ^ + i2 5? 


Id,. , . du 7 

■» ■ s~ s; (hl%) + hl sr 




where = u^a^y), Up = Up(a,0,7), and u y = u 7 (a,0,7) are 
the components of the displacement vector of the poinlj (a,0,y) on 
the axes of the orthogonal trihedron, the vertex of which is at the 
point (a, 0 , 7 ) and the faces of which coincide with the planes 
tangent to the surfaces a = constant, p * constant, and 
7 = constant. The positive direction for the displacements cor- 
responds to the direction of increase of the coordinates a, 0, 7. 


Equations ( 2 . 1 ) are oDtained from the general formulas of the 
theory of elasticity given for example in the hook of L. S. Leiben- 
son (reference 1) for 113 = 1. 

The equations of equilibrium of the general problem of the the- 
ory of elasticity in the coordinates of the shell are presented in 
the form 
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( 2 . 2 ) 


where O and T are the normal and tangential stresses . The sub- 
scripts denote, for normal stresses, the direction of the outward 
(in the direction of increase of the corresponding coordinates) 
normal to the corresponding surface; for the tangents, they denote 
the surface of action of these stresses taken in pairs from the 
conditions of their reciprocity. The components of the stress 
tensor are considered positive if, when applied to the surface with 
positive outer normal, they are directed toward increasing coordi- 
nates. The magnitudes p a , Pp, and p^, in equations (2.2) axe 
the components of the vector of intensity of the volume forces. 
Equations (2.2) are obtained from the general equations of the the- 
ory of elasticity given for example by Love (reference 2) for 
h 3 = 1. 

In the theory of shells, the stresses O a , Op, r a p * Tp a , 
applied normal to the section and lying in a plane tangent to 
v = constant, are determined from the six stress components 
expressed in terms of the strains. The remaining three components 
of the stresses axe found from the conditions of equilibrium. 

From Hooke's iaw, only the three relations referring to the 
stresses C^, Gp , and T a p are retained; these relations are 
given in the form 

(? a = (k+2|i) A - 2\i (epp+e^) 

0p = (X+2p) A - 2|i (e^+e^) 

T aP = T j3a = l-^ap 



(2.3) 
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where A is the volumetric dilation and \ and p are the coef- 
ficients of elasticity of Lame. 

The equations of equilibrium (2.2), on the basis of equa- 
tions (2.3) and (2.1) after a number of transformations using, where 
required, the relations (1.3), (1.4) and (1.5), are given in the 
form 


<K * 2|l) EJ ST 2 ** H7 + 2 “ ABK “a sr) + h i 

. . . 1 d& 1 dX d / 1 

(X+2u) hi +2u ^ + 2uABKu p - 2p 3^ )+ h 2 


- 2 ( X+2p ) (H+Ky) ABA + 2 p 


d_ 

Ba 


( Bk 2%) + ^ 


(Ak]_Ug ) 



4uAB(H+Kr) 



0 


(2.4) 

where K = K(a,f3) and H = H(a,p) are the Gaussian and mean cur- 
vatures of the coordinate surface, respectively. 

K = k x k 2 
H = g- ( kq+k 2 ) 



The volumetric expansion and the normal component (the pro- 
jection on the normal to the surface y - constant) of the ele- 
mentary rotation of the shell are denoted by A = A (a, (3,?) and 
2X = 2X(a,0,y), respectively. In the following discussion 2X 
shall be denoted simply the normal rotation. The volumetric expan- 
sion and normal rotation are determined in terms of the displace- 
ments u (X , Up, and u^ in the coordinates of the shell by the 

formulas 
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Equations (2.4) in orthogonal coordinates for h.3 = 1 are the 
general equations of equilibrium of an elastic body. The first two 
of these equations express the tangential equilibrium of the three- 
dimensional element dadpdy/hphg of the shell, that is, the equi- 
librium of this element in the plane tangent to the surface 
7 = constant. The last equation refers to the equilibrium of this 
element in the direction of the outer normal to the surface 
7 = constant. Equations (2.4) differ from the equation of the gen- 
eral problem of the theory of elasticity in displacements or 
strains in the fact that each of them contains both static and 
kinematic magnitudes. 

5. Displacements and Strains of the Shells . - The theory of 
shells is based, as is known, on the hypothesis of Kirchhof f -Love 
according to which a rectilinear element normal to the middle sur- 
face of the shell remains, after deformation, rectilinear normal to 
this surface and of the same length. This hypothesis is equivalent 
to the assumption 

e^y “ ~ ^77 ^ ^ (3.1) 

and leads, for the displacements Uq,, up, and u^ of an arbi- 
trary point (a, 0,7) to the formulas 

% - (MV> « - 5 

Up = (l+k 2 r) V - I | 

Uy = W 

where u = u(a,0) and v = v(a,0) are the tangential displacements 
(in the direction of the tangents to the lines 0 = constant and 
a = constant) of the point (a,0) of the coordinate surface, and 
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w = w(a,|3) is the displacement of the same point in the direction 
normal to this surface, positive in the direction of increasing 
coordinate y (fig. 2) . 

For the purpose of presenting a more accurate theory valid 
not only for shells of medium thickness hut also for thick shells, 
another hypothesis, which is a generalization of that of Kirchhof f - 
Love, will he used. 


It shall he considered that each of the three components Uq,, 


Uj^ y ClIJ.1 

setting 


and u v is represented as a function of 


hy a linear law. 


(U, 1? ) U.I| 


= (l+k 27 ) v - I 


ly = W + 7W* 


(3.3) 


where u, v, and w have the same values as in equations (3.2); 
w* = w*(a,p) is a magnitude that depends, like u, v, and w, 
only on the two variables a, (3 and is the relative elongation of 
a normal element of the shell (constant under the assumption made 
here over the entire length of this element). It is easy to see 
that with equations (3.3) equations (3.1), which express the funda- 
mental hypothesis of the present theory of shells, do not apply. 
With the introduction of the deformation of elongation w* of a 
normal element of the shell, all the six components of the strain 
tensor (2.1) receive values different from zero. 

Equations (3.2) establish the kinematic model of the deformed 
state of the shell. This state, in the general case, is made up of 
two states of which the first is determined only on the tangential 
displacement u,v of the point of the coordinate surface (w, w* 
in this case being equal to zero) and in the second only hy the 
normal displacement w and the elongation w* (u,v in this case 
being equal to zero) . The deformation of the shell determined only 
hy the tangential displacements u,v shall he denoted the tangen- 
tial deformation for briefness. This deformation is characterized 
hy the fact that an arbitrary point of the surface y = constant 
after deformation does not go beyond the limits of this surface as 
a two-dimensional space. An elementary layer of the shell d y for 



8 


RACA TM 1241 


a tangential deformation does not change its shape and position in 
space and undergoes deformations of length and shear at the sur- 
face 7 = constant as a two-dimensional space (in the general case 
for K / 0 non-Euclidean) . A deformation of the second kind 
determined only by the normal displacements w and the elongation 
v* will be called a normal deformation of the shell. For this 
deformation, an arbitrary point (a, (3, 7) of the surface 7 * con- 
stant passes with respect to this surface into the third dimen- 
sion. A normal deformation is accompanied by a change in shape 
of the surface. 

In setting up the kinematic model determined by equations (3.3) 
for all six components of the deformation tensor, by virtue of 
equations (2.1) and (1.2), a definite law of variation with thick- 
ness of the shell is obtained. 

Representation of e^, epp, and e a p in the form of series 
in the variable 7 gives 

e oa = e l + ^ X ln y" 

e 3P " e 2 + X X 2n ^ 
e aP = a,+ X Tn 7n 

(n = 1,2,3, . . .) 

where the coefficients of the series 6j_, Xq n , X 2n , 

and v n each depends only on the displacements u, v, and v* 
of the point (a,p) of the coordinate surface. By substituting 
the displacements %, up, and u^, determined by equations (3.3) 
on the right-hand sides of the corresponding equations (2.1) and 
then by representing the magnitudes hq, h 2 , and their ratios 
(direct and inverse) in the form of series 
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1 /t i ,22 ,3 3 \ 

hi = i (1 - k]_7 4- k]_ 7 - k]_ 7 + . . . ) 


, 1/ni ,22 ,33 x 

&2 = ~ (1 - k 2 7 + k 2 7 - *2 7 + • • • ) 

.D 


~ “ | [} - (k!-k 2 )(7 - k l7 2 + k x 2 7 3 ' • • • j) 

2 

J“ = | |l + (k x - k 2 )(r - k 2 r 2 + k 2 2 7 3 - ... 2 I 


(3.5) 


and referring to the last two of relations (1.3) and relations (1.4), 
after a number of transformations for the coefficients of the ser- 
ies (3.4) the following equations are obtained; 

€q = i^]i+-i_^V4- kqw 

A 3a AB 3(3 

e o = _!_^2.u + i^Z + k 2 w 
AB 3a B 3(3 

00 = A JL/hA + B 3_/u\ 

3 b$\A J A 3a yB j 

Xln = (-l) 11 " 1 ^" 1 ^ E + ^ I-ki 2 w-i JL/l - JL M^k-kW*] 

_Sa A Sp B A Sa \^A Sa J ^g2 S 3 S 3 


X 2n = (_1) k : 


1 n-l[^k 2 u 3 k 2 y 2 1 3 / 1 3w\ 1 3 b 3w 


2 |_3a A + 3(3 B 2 W B 33 \b 3(3 J 3a 3a 

c~ 

%- (-l) n_1 iCkl-ks) fks 11 - 1 1 3 M - tL 11 - 1 5 S M 

L L B °P \A/ A 3a \ b ) 


+ k 2 w* 


L n-1 B _ 3 _ / v 
A 3a l B 


, n-1 . n-1 

kq + k 2 

3 2 v 

1 3b 3w 

1 3a 3w ^ l 

AB 

\Sa S3 

A 3a 3(3 

B 3p 3a / j 


(3.6) 



10 


NACA TM 1241 



In the following discussion, formulas for the volume expansion & 
and the elementary rotation 2X will he required . When these mag- 
nitudes are also represented in the form of series in powers of 7 


4 = 4 o V" 

*-*(>♦£ X„7 n 


(3.8) 


Then by making use of equations (2.6) and (3.5) after a number of 
transformations, using relations (1.3) aud (1.4) for the coefficients 
of the series (3.8), the following formulas are obtained: 
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A 0 = €q + € 2 + V* 


AB 


(Bu) + A. (At) 
oa ^3 


+ (k^+kg) w + w* 


A„ = x 


n = ■‘•In 


x 2n - (-D"- 1 (A' 1 ^ A 1 

da oa / A 


X kl 

n-1 


n-1 Sk l 


+ k 2 


n-1 Sk 2 \ v 


Sp T “ z & J i ' 111 

1 5 /I 5w\ t 1 5 b ow 
B 5|3 \B 5p J 5a 5a 


n-1 


1 _5_ 

f 1 5w\ 

A 5a 

\A 5a y 


1 5a 5w 


(ki n+1 + k 2 n+1 )w + (ki n + k2 n )w* 


X 0 = 


2AB 


A (By) - A (Au) 

5a 5p 


X n = (-1) 


k, -ko r 


n j -nq-a-g 


n-1 a 5 


- \ + k- 


L 2 B Sp l A 


n-1 b 5 
L A 5a 


(S' 


n-1 


- k. 


n-1 


2AB 


2 ( 5 w l 5a 5w 1 5b 5w 

5a 5p . A 5p 5a B 5a 5p y 






(3- 


J 


4. Analysis of Kinematic Relations. Corrections and Additions 
to Theory of Love . - Equations (5.4) and (5.6) for the components 
of the deformations have a common character and were obtained in 
correspondence with the hypothesis (3.3) assumed for the displace- 
ments. For w* = 0 from equations (3.3, 3.4, and 3.6), there are 
obtained equations for the displacements and the deformations of 
the shell having an inextensible normal element and following the 
hypothesis of Love (3.2). The magnitudes Xqq, Xgq, and Tq, 
defined by equations (3.6) for n = 1 and w* = 0 and the first 
two representing bending deformations (variation of the principal 
curvatures kq and kg) and the third the torsional deformation, 
differ from the corresponding magnitudes Xq, X 2 , and T , which 
were used by Love. By setting w* = 0 in the last three equa- 
tions of (3.6), 
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Sa A Sp 


1 I - k x 2 v - 

A Ba \A Ba 


B 


ok. 

* 21 ’ £ 
T 1 = (ki-kg) 


Bk c 


2 u UA 2 t ,2 1 3 

— + x- k, w • - 

A Bp 


B 


B Bp 


1 Ba By 

AB 2 dp dp 


/ 1 Bv\ 1 

\B Bp J A 2g 


BB By 
Ba Ba 




[(4 J) 


Ai(! 

i \_ B _B_ /v\ 

__L 

d 2 w 

1 Bb By 

_ 1 Ba 

B Bp Vj 

aJ A Ba w 

AB 

^Ba Bp 

A Ba Bp 

” B Bp 




From these equations, it follows that for tangential deforma- 
tions (in the case w = 0) the changes in curvature X^, and 
are determined as linear algebraic expressions relative to the dis- 
placements u and v with coefficients proportional to the partial 
derivatives of the principal curvatures k^, and k2 of the unde- 
formed surface. The expression for the torsional deformation will, 
as is to he expected, he symmetrical with respect to the coordinates. 
The same properties, as seen from equations (3.6), are possessed 
also by the remaining components X^ n , Xg n , and T n for 
n = 2, 3, 4, . . . . In particular, for the spherical shell a 

result is obtained that generalizes in a certain sense the theory of 
the bending deformation of a plate as based on the hypothesis of 
Kirchhoff. This result can be formulated in the following theorems; 

Theorem I . - The deformations of elongation and shear e<xa, 
epp, and e&p and the volume deformation A of a spherical shell 
in the case of tangential deformations (that is, for w = 0) are 
uniformly distributed over the thickness of the shell (do not depend 
on y) and are determined only by the deformations of elongation 
and shear e €g, and u) of the middle surface. An exception 
to the uniform distribution of the magnitudes e^, e pp> and e a p 

over the thickness of the shell arises only as a result of normal 
displacements . A change in the shape of the spherical shell char- 
acterized by the parameters of the change in curvature X-q, Xpq, 
and t-j_ is due only to the normal displacement w. 

Theorem II . - The normal rotation 2X of the spherical shell 
is determined only by the tangential deformation (the variables u 
and v) and remains constant over the thickness of the shell. In 
the case of normal deformation, the normal rotation 2X is equal 
to zero. 

This result, obtained on the basis of the analysis of the gen- 
eral formulas of the preceding section for the spherical shell, may 
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be arrived at directly in the following manner. It is assumed that 
the spherical shell as a deformed body is, at one of its bounding 
surfaces (for example, the inner), in contact with a rigid spheri- 
cal base so that an arbitrary point of the shell can be freely dis- 
placed along the surface of this base without going outside the 
limits of this surface. Such a model corresponds to the case of 
tangential deformation of the shell. Now at the point (a, (3 ) 
some normal section of the shell, in general arbitrary with respect 
to the chosen coordinate lines a and (3, is assumed. The lin- 
ear normal element, as the shell passes into the deformed state, 
remains, by the Love hypothesis, normal to the base surface and 
takes on a new position determined by the rotation of this element 
with respect to the center of curvature (in the case of a sphere, 
common for all normal elements). Let M]_ ' Mg ' be the projection 
on the plane of the chosen section of the element M-jMg after 
deformation (fig. 3). Further let £ = M]M^' denote the projection 
on the plane of this section of the vector of the total displace- 
ment of the lower point of the element. Then the displace- 

ment i y = MM' of an arbitrary point M of the element M-jMg in 

the plane of the chosen section will be equal to 

e 7 = i (1 + kr) (4.2) 

where k = l/R is the curvature of the inner surface of the shell 
and y is the distance of this surface to the point M considered. 
The corresponding elongation of the tangential element 
ds = (l + k/) Rd9 is determined by the equation 

6 = 5s^ = R(1 + k7) h P 1+k7) *] = (4 ’ 3) 

From equations (4.2) and (4.3), it follows that whereas the 
tangential displacement £ 7 of the spherical shell is a linear 
function of the coordinate 7 , the deformation of elongation e 
does not depend on 7 . The same result can also be obtained 
directly from equations ( 2 . 1 ) for the deformations e^, 
and e a p for the values entering the following formulas : 

ki = kg * k = constant 

*1 - ■ rrw 

u a = (1 + k7)u 

^ = (l + k 7 )v 
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In the same way on the basis of the second of equations (2.6), 
the second theorem can be proven. 


In constructing this theorem, Love represents the components 
of the deformation e^, e^ , and e a p in the form of linear 
expressions relative to the parameter y. 


e aa = e l + X l? 
e PP = £ 2 + X 2? 
e ap = “> + J 


(4.5) 


and for the parameters of the change in curvature Xq, X2, and T 
gives the equations 


v 1 d /, 1 1 dA 

Xl ■ I £ (k i u) + 


AB dp 


(k 2 v) - ~ J 7 {- i“ 
c A da \A da 


) 


1 dA dw 
AB 2 dp dp 


_1_ dB 
AB da 


(kqu) + - ~ (k 2 v) — 

1 B dp 5 6 B dp \B dp/ a 2 : 


T ■ I £ ( ^ T) 


1 dv ^ 1 d / 1 dw 

A da A da \B dp 


1 dA dw 
A 2 B Sa 




dB dw 
B da da- 


r ( 4 -6) 


J 


These equations differ essentially from equations (4.1). The 
magnitudes Xq, X 2 , and T determined by Love as coefficients of 
the second members of equations (4.3) are in contradiction to the 
theorems just proven for the spherical shell. The difference noted 
here in the determination of the magnitudes Xq, X 2 , and t by 
equations (4.1) and (4.6) is explained by the fact that Love and 
other authors (in particular, Timoshenko (reference 3)), following 
Rayleigh, start from the assumption of the inextensibility of the 
middle surface. This assumption stands in certain contradiction 
with the geometry of extensible and flexible surfaces. 

In recent years a number of papers have appeared that refine 
to a greater or less extent the theory of thin shells of Love . 

The most interesting and original of these are the investigations 
of Krauss (reference 4), N. A. Kilchevsky (reference 5), and 
A. I. Lurie (reference 6 ). 

5. Fundamental Differential Equations of Equilibrium of Elastic 

Shells . - The general equations of a shell possessing deformable 
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normal elements is obtained by reducing the three-dimensional prob- 
lem of the theory of elasticity to two-dimensional starting from 
equations (2.4) , retaining in the series (3.8) the first three 
terms, and applying the principle of Lagrange corresponding to the 
kinematic hypothesis (3.3). 

An element of the shell AB5dad(3, having an infinitely small 
area ABdad(3 on the middle surface and a finite length 5 equal 
to the thickness of the shell, possesses according to the kinematic 
model seven degrees of freedom; namely, six degrees with respect to 
the displacements of the element (three linear and three angular) 
in space as a rigid body, and one degree characterized by the 
change in length of the element. Corresponding to these degrees 
of freedom seven equations of equilibrium must be obtained. Of 
these equations, the first six refer to the equilibrium of the 
element in space as a rigid body and the seventh may be obtained 
by equating to zero the work of all the external and internal 
forces of the element AB&dadp against displacements and deforma- 
tions corresponding to the unit elongation w* = 1. It should be 
noted that the equations of equilibrium of an element may also be 
obtained on the basis of the principle of virtual displacements by 
equating to zero the sum of the work of all the forces (in the 
given case only the external, because the element is considered as 
a rigid body) for each of the six possible unit displacements. 

By the method assumed here, one of the conditions of equilib- 
rium of the element as a rigid body, namely the condition corres- 
ponding to the rotation of the element about the normal to the mid- 
dle surface and given in the theory of Love, the sixth nondiffer- 
ential (relative to the shearing forces and torsional moments) is 
satisfied identically because of the relation T a p = Tp a used in 
deriving the general equations (2.4). 

Thus, starting from equations (2.4) and applying the principle 
of virtual displacements, it will be necessary to obtain for an 
element of the shell only six equations, one of which (called above 
the seventh) according to its physical meaning represents the gen- 
eralized condition of equilibrium of the element ABSdadp having 
a strain w* expressed as a function of y. 

Substituting in the left sides of equations (2.4) the dis- 
placements Ug., Up, and u^, according to equations (3.3), the 

volume dilation and the normal increment A and X according to 
equations (3.8) (in which it is necessary to retain only the first 
three terms, that is, to y ^ inclusive and reject the others) and 
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hi and h2 according to equations (1.2), three expressions are 
obtained each of which contains terms with powers of y up to the 
third inclusive. 


The three magnitudes thus obtained represent, according to 
their physical meaning, the components along the axes of the mova- 
ble trihedron on the surface 7 = constant of the vector of the 
external force acting on the three-dimensional element of the 

shell d y and expressed in terms of the kinematic magnitudes 

h l h 2 

v, w, w*, A q , A^, & 2 > Xq, *i> and *2 8X1(1 ^atic magni- 

tudes T^, t p, Cy, V a > Pg> 311(1 Py In passing to the two- 
dimensional element of the shell ABSdadp, the work of all the 
forces acting on this element and determined in this manner must 
be equal to zero on each of the five possible displacements as a 
rigid body. 


Corresponding to these displacements and by virtue of hypothe- 
sis (3.3), each of the first two equations of (2.4) must be by 
d y and ydy and the third by d y, integrated with respect to y 

between the limits 7 = - i 6 to r e + i 6, and the result 

2 2 * 


equated in each case to zero. Thus five equations are obtained 
containing in addition to terras with the kinematic magnitudes u, 
v, w, w*, Aq, Aq, A2, Xq, Xi, and X2 also terms with the 
transverse forces N-j_ and arising from the tangential 

stresses T and t^. 

rl 8 1 


Ni - «• 


No 


\J 


yq 

1 h 2 

’I 5 


dy 


1 

A 


\J 


n + I 5 

IlP d.7 
i. h i 
? 8 


J 


(5.1) 


In order to obtain the sixth equation corresponding to the lin 
ear strain of w* of the normal element of the shell, the left 
side of the third equation in (2.4) must be multiplied by ydy 
and the integral of this expression between the limits 

7 « - i 5 to y = + 8 equated to zero. When it is remembered 

that 
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n 


v - 


^(^) rd7 


il l h 2 


Is 


1 R 

2 S 


1 8 


sj 


1 R 

2 6 


h l h 2 


dr (5.2) 


where the first term refers to the work of the external and the 
second to the work of the internal normal forces of the element AB5 
for the normal displacements u^, = yw* for w* = 1, determining 
Typ , and Oy in terms of the deformations by the equations 


v - ^ V s ) \ £■ 

- ^-^ 2 > I If 


>1 

> 


Oy = X(A 0 +A 1 r+A 2 r 2 ) + 2|aw* J 


(5.3) 


and representing the remaining terms of the third equation of (2.4) 
in the form of a finite series in powers of y, an equation is 
obtained in which the unknown will be only one of the kinematic 
magnitudes . 

Thus there are six equations with respect to 12 functions, 
the four basic functions u, v, w, and w*, the six func- 
tions Aq, Aq, A 2 , Xq, Xq, and X 2 giving the coefficients 

of the first three terms of the series (3.8) and the two transverse 
forces Nq and N 2 . 

These equations, upon eliminating the forces Nq and Ng, 
reduce to four equations. Neglecting the small terms with k^-i6 2 /l2, 
k / 12 y and k-jkgB /12 in the expressions 1 + k ]_5 / 12, 

1 + k 2 2 5 2 /l2, and 1 + kqk 2 S 2 /l2 finally 
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+ dp \B 1 dp / 
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d 
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4 o "a - rv' I . 1 T NL' / 4" ~v ' j -x r '"' I 


12 AB[3a \^a Sa / Sp yB dp y 


- 2{JIV*+| Z* * 0 


(5.4) 


where X = X(a,p), Y = Y(a,p), and Z = Z(a,p) are the components 
on the axes of the movable trihedron of the vector of surface 
intensity of the load computed for the stresses t a y , T^y , and Oy 

on the boundary surfaces 7*^8, 7 = - 5 and for given volume 

forces P a , Pp, and p^, by the formulas 
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(5.5) 


The magnitudes = n^(a^p) and m^ = n^(a,p) are the 
moments of all the forces (surface and volume) relative to the 
axes a, (3 of the movable trihedron of the middle surface: 



Finally the magnitude Z* = Z*(a,f3) is the new generalized 
static magnitude corresponding to the elongation of the normal ele- 
ment and determined by the equation 

+ 5 + \ & 

p y ° y y 

— — yd 7 + — — 

1 * h l h 2 h l h 2 1 c 

2 ° “ 2 ° 
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In the case where only surface forces act on the shell, the 
first components in equations (5.5) to (5.7) drop out. 

To the equilibrium equations (5.4) must he added the equations 
for the components Aq, A-^, and A 2 of the volume dilation and 

X Q , and Xg of the normal rotation. According to equa- 

tions (3.9), 


At =2 


A n = — (Bu) + ( Av)"l 

o AB Loa dp J 

dB u dB v\ 1 | d / B dw\ _d_ ( A dw\ 
v da A + dp B/” AB da\A da/ dp \B dp/ 

2 [da 


+ 2Hw + w* 


- (k^l+h^2)'w + 2 Hw* 


/, 2 , 2 , U d/,2 , 2 \ V 

(k x +k 2 ) A + d? k 1+k 2 B 


v i. _d_ / 1 dw d 1 dA dw 

•*- _A da V A da / £)p dp 


w 

3_ 


k. 


!L _d_ / 1 dwN 1 dB dw~ | 
: B dp \B dp / j^B da daj 


+ (k^+k^ 2 )w - (k^^+k 2 2 )w* 


X n = 


0 ~ 2AB da 


S (Bv) - A- (Au)l 




J 


- L 


B _d_ 

A 5a 


A + ± k (A 

B ) B dp \AJ 


X 2 = L 


k, B A (l\ + kp A i_ ( 

A + S 2 w 

1 dA dw 

1 dB dw 

1 A da \BJ B dp v 

\J dadp 

A dp da 

B da dp 


(5.8) 


In equations (5.4) and (5.8), 


K = k x k 2 
H = l (k 1+ k 2 ) 

L - | ( k l~k 2 ) 




(5.9) 


J 


Equations (5.4) together with equations (5.8) form a complete 
system of differential equations of the shell. To these equations 
the Boundary conditions for each particular case must he added. 
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For this purpose , the internal generalized forces must he 
determined. These forces, corresponding to the degrees of freedom 
of the normal deformable element, will in any normal sections of 
the shell consist of the tangential (normal and shearing) forces T 
and S acting in the plane tangent to 7 = 0 and corresponding 
to the displacements of the element parallel to this plane, the 
transverse force IT directed along the normal to 7 = 0 and cor- 
responding to the displacement of the element along the normal to 
the middle surface 7=0, the bending and torsional moments G 
and H corresponding to the angular displacements of the element 
with respect to the tangent axes of the movable trihedron and, 
finally, the new generalized (statically equivalent to zero) trans- 
verse force N* corresponding to the elongation of the normal 
element . 

All of these forces, with the exception of the transverse 
force K, can be expressed in terras of the fundamental kinematic 
magnitudes u, v, w, and w* by setting up the work of all 
(tangential and normal) stresses of the normal section considered 
over unit displacements of the normal element, translational in 
the tangent plane, translational in the direction of the normal 
to the middle surface, angular relative to the axes in the tangent 
plane, and in the displacement of the points of the element u^, = w*7 
for w* = 1. 

For the internal forces on the two basic normal sections 
a = constant and (3 = constant (fig. 4), the following equations 
are obtained; 
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In these equations, 


1 5u l 5 a , 

£-| = + V +• k]W 

A 5a AB 5b 
1 SB 1 5v . 
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1 u og l v , c 
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/5 2 w 

l 5 b 5w 

l 5 a 5w\ 

U) 

” A 5a (.By 

" AB 

\5a5p 

" A 5a 5 b 

’ B 5b 5a/ 


(k -k ) k A _5_ /u\ _ k B 5 / vN + k l +k 2 / 5 2 y _ 1 5b 5w 1 5a 5w \ 
1 2 |_ 2 B 5p \AJ 1 A 5a \bJJ AB \ 5a5p A 5a 5p B 5p 5a/ 


(5.11) 


Depending on the character of the problem, the boundary con- 
ditions may be purely kinematic, purely static, or of the mixed type. 


In the case of kinematic conditions for the normal element of 
the shell and the boundary surface, there must be given in the boun- 
dary surface three displacements of the midpoint of the element 
along three mutually perpendicular directions, the angle of rotation 
of the element relative to the tangent to the contour curve of the 
middle surface and finally the normal displacement of any other 
point of the element. Altogether there will be five independent 
kinematic conditions, which together with the fundamental equa- 
tions (5.2) and (5.8) make the problem entirely determinate. 

If the boundary conditions are given in terms of stresses, 
there will be in this case five independent conditions, the four 
usual conditions of the moment theory and one with regard to the 
generalized (statically equivalent to zero) transverse force N*. 
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6. Shells of Medium Thickness . - By neglecting in eq.ua- 
tions (5.4) and (5.8) the small terms with tangential displace- 
ments u and v, which contain as factors the products of the cur- 
vatures k]_ and kg , and their derivatives for the shell of medium 
thickness 


w* 

\ = 

M- * 


0 


Ey 

1-u^ 
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2(1+ U) 


'N 




J 


(6.1) 


and neglecting, in correspondence with this the last equation of (5.4) 
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(6.3) 


For the internal tangential forces and the moments, 


Tl = 1^2 f 1 +U£2 " 12 ( k l' k 2) X lJ 

T 2 = [e 2 +ue ! + (k 1 -k 2 )X 2 J 

S x = -M Fo>+ §0- (k 2 T 1+ T 2 )l 

1 2(l+u) L 12 2 1 2 1 
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where the strains £]_, £ 2 , w l X-^, T i> T 2 are deter- 

mined by equations (5.11), in wiiich (fourth and fifth) w* must 
be set equal to zero. 


Equations (6.2) to (6.4) refer to shells of medium thickness 
and were obtained with an accuracy up to terms with S 3 /l2 in 
strict correspondence with the fundamental assumption (3.2) because, 
in the first place, the given equations (4.1) for the components 
of the changes in curvature X-^, X 2 , and Tj_, in contrast to 

equations (4.6) of Love are accurate; and in the second place, these 
equations contain additional terms arising from the moments, namely, 
the terms with A 2 S 3 /l2, and X 2 S 3 /l2, which are absent in the 
present moment theory constructed on the basis of the magnitudes e^^, 
epp, and e a p in the form of equations (4.5) and not in the form 


e oa “ £ 1 + *117 + X 12X 


= £ 2 + X 21 7 + X 22 7 


a(3 



J 


(6.5) 


which lie at the basis of the theory given here. This theory and 
the more general one given in the preceding section and referring 
to thicker shells is in full agreement with the fundamental theorems 
of the theory of elasticity, in particular, with the theorem of 
reciprocal work, which, as shown later for the example of a cylin- 
drical shell, does not correspond to the theory of Love. 


7 . General Technical Theory of Thin Shells . Two methods of 
Solution of the Problem. Generalization of the Maxwell and Sophie 
Germain- Lagrange Equations . - For thin shells, further simplifi- 
cation of equations (6.2) to (6.4) is possible. Eliminating from 
equations (6.2) the functions Aq, A^, A 2 , Xq > X l> and X 2 

with the aid of relations (6.3), three equations are obtained in 
the three functions u, v, and w. The first two of these equa- 
tions, when multiplied by 5, will each consist of terms propor- 
tional to the thickness 5 and terms proportional to magnitudes 
consisting of the product of 5 3 /l2 by the curvatures, k-^ and k 2 
or the -derivatives of these curvatures . 


In the third equation, in addition to the terms of this type, 
there will enter a term proportional to B 3 /l2 and independent of 
the curvature of the shell. 
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The rather extensive theoretical and experimental investiga- 
tions made by the author show that for thin shells the relative 

,3 *3 

k 2> 
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12 Sa 12 Sa 12 dp ' 12 dp 

mental equations are factors of 3econd-order values for the dis- 
placements u and v. Without sensible error these terms, as 
shown in the work on cylindrical shells (reference 7) and thin- 
walled rods (reference 8), can be neglected. Correspondingly, 
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(7.1) 


where because of the assumptions made, 
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After certain transformations and simplifications of equa- 
tions (6.2), 



NACA TM 1241 


29 


2 

V e 2 V e 2 9+ 2V e 2 (Hw) - (l-p)(HV e 2 -LV h 2 )v - - ^ ± [|j- (BX) + (AY] 

-a* e V (I-vKHV^-LV, 2 ),. (1-U) £[£ (* 2 |)-| (h g)] - 
i' ? aV''-2[2H 2 -(l- l .)K]v=-^Z 

(7.4) ' 

vhere cp, \p, arid v are the required functions of the displace- 
ments and are invariant (relative to the directions of the coordi- 
nate curves a and p at a given point of the surface) magnitudes; 
V e 2 and V h are the differential operators of the second order 
of the elliptic and hyperbolic type: 

v 2 . i. [i- fz £-) + (± l-'ll , 1 
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A da \A da/ B dp \B dp / ^2-g da da, ^2 dp dp 


The mixed operator HV 2 - LV 2 is defined by the equation 

e n 
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(7.6) 


Differential equations (7.4) form a complete system of equa- 
tions in the three fundamental displacement functions q>, f, 
and w. 


These functions according to equations (7.3) determine the 
vector of total displacement of the point (a, (3 ) of the middle 
surface and therefore hy virtue of equations (7.2) and (7.1) all 
the deformations and the internal axial forces and moments of the 
shell. 

Equations (7.4) are thus the fundamental equations of the the- 
ory given here for thin shells and permit solution of the problem 
of the equilibrium of elastic shells of small curvature by the 
method of displacements . 


The theory of thin shells can also be presented in another 
more compact form, namely, in the form earlier proposed of the 
mixed method by introducing only two functions, the stress func- 
tion $ and the displacement function w. Setting (for X = Y = 0) 


T _ 1 d /I <&) 1 dB d$ 
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1 dA d$ 
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(7.7) 
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and bearing in mind the analogous equations (7.2) for X^, Xg, 
and ir and equations (7.1) for G^, G g , and = - Hg, the 
general equations of equilibrium and deformations of the moment 
theory of shells are represented in the form of two symmetrically 
constructed differential equations : 

h V eV* - » 7 e 2 -Wh 2 )’' - ° 

- - E65 %\ Z '* * 2 
12(l-u ) 



These equations are a generalization of the equation of Max- 
well for the two-dimensional stress state of a plate and the equa- 
tion of Sophie Germain - Lagrange for the case of the bending of a 
plate, inasmuch as for k^ « kg * 0 (the case of a flat plate) 
they break down into the well-known equations from the theory of 
elasticity 



v e 2 $ 


0 


7 2 V 2 v 
© © 


12(1-1^) 


ES V 


(7.9) 


O 

in arbitrary (for V Q determined by equations (7.5)) orthogonal 
coordinates . 


If in the second of equations (7.8) the term with 8^/12 is 
neglected, the fundamental equation of the momentless theory of 
shells results : 


(HV e 2 -LV h 2 )$ = Z 


(7.10) 


After determining the functions $ and w, the forces 
and T£ are found from equations (7.7), the moments G^, Gg, 
and 1 from equations (7.2) and (7.1). 

These forces and moments will satisfy the equations of equil- 
ibrium 



32 


NACA TM 1241 


C < OT i> - I - 1 <**> + s i | + “A + "* ■ 0 

(AT P ) - T, ^ + A (BS, ) - S P ^ + ABk 0 N 0 + ABY = 0 
dp * 1 OP da 1 * da 




2 2 


- (k-jTj+kgl^) + A 


d (BN X ) + A (AN 2 ) 


da 

d 


+ Z = 0 


A (BH-, ) - H„ A - A (AG p ) + Gn ^ - ABNo = 0 
da 1 2 da dp 2 1 dp 2 

A (AHo) - H-, ^ + A (BGi) - Go— + ABNi = 0 
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+ S 2 + kj^HjL + k 2 H 2 = 0 


(7.11) 


for X = Y = 0 with an accuracy up to the terms ABk-jN^ and 
ABk 2 N 2 in the first two equations and the term kjH^ + k 2 H 2 in 
the last equation, which as magnitudes proportional to k^8 3 /l2 
and k 2 8 3 /l2 (hy virtue of the fourth and fifth equations of (7.11) 
and the relations (7.1) and (7.2) for the moments) are taken equal 
to zero . 


In neglecting in the first two equations the terms with k-jN^ 
and k 2 N 2 , an error is admitted of the same order as that in the 
general theory in replacing the last of equations (7.11) hy the 
approximate relation * - S 2 . 


For the transverse forces N]_ and N 2 , 
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(7.12) 


which constitute a generalization of the well-known equations of 
the theory of the bending of a plate. 
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8. Circular Cylindrical Shell. Particular Cases . - For the 
coordinates of a and p, the distance to the point considered 
along the generator and the transverse arc, respectively, of the 
middle surface are taken. 


ko « 


Then, evidently A = B = 1. Equations (5.8) for k-i = 0 and 
1/B « constant assume the form 
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2 \£iadp 2 dp/ 


> ( 8 . 1 ) 


J 


Equations (5.4) may be presented with the aid of table 1, con- 
taining the differential operators. 
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In this table X, Y, Z, and Z* denote free terms depend- 
ing on the internal forces ( 1 % and mp are taken equal to zero ) . 
The differential operators referring to the corresponding func- 
tions are indicated. These operators form a symmetrical differen- 
tial matrix. The elements of the matrix symmetrical with respect 
to the diagonal terms hare the same express ions , a consequence of 
the theorem of reciprocal work. 


For k 2 = 0, equations (8.2) break down into the following 
equations (likewise symmetrically constructed): 

(\ 0 \ d 2 u c) 2 u ^ S 2 v % dw* 1 „ „ ] 

(A +2(i ) — — + |j. — — + ( A+p ) — N + A - — + — X - 0 




d 2 v ow* 1 v A 

— s- +A - — + - Y = 0 

;w 2 Sp 6 


A ^ + A — - n? 2 w* + (A+2(i)w* - — Z * 0 


and the equation for the bending of a plate 


V 2 V 2 w = 


(A +2(1)5" 


Equations (8.3) are, in a certain sense, a generalization of 
the problem of the two-dimensional stress state of a plate and per- 
mit determining the stresses and the deformations of the plate 
under the action of two mutually balancing concentrated forces 

applied on the planes y = + i 5 and y = - 5 and acting normal 

to the middle surface. In the case of the homogeneous problem, 
the first two equations of (8.3) may be satisfied by setting 

u = - 

oa 

y - -A V 2 $ 

ap 

w* = (A +2(i ) V 2 V 2 $ 

where $ =$(a,p) is an arbitrary function. The last equation then 
becomes 
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y 2 y 2 y 2 $ , V 2 V 2 $ m Q (8,6) 

\+2p 

and therefore the function - 4 (X-Hjl)/(\+2|j.) is hiharmonic . 

The magnitudes u, v, and v* determine the strains e aXL } 
e pp , e^ , and e^ and therefore the stresses O a , Op , and 
T a p . The remaining stresses, sis in the general case of the shell, 
must he found from the condition of equilibrium. 

In the same manner as the peirticular case of equations (8.2), 
there can he obtained the fundamental equations for the circulsir 
arch with account taken of the extensions! deformations of the eirch 
in the direction of the normal to its axis. In this case, the dis- 
placement v must he considered equal to zero and the remaining 
magnitudes u, w, and w* considered only as functions of p. 
Equations, which generalize the well-known equations of Boussinesq, 
are obtained. 

For 


w* = 0 



** 2(1+1?) V 

(where 1? is the Poisson coefficient), equations for the circular 
cylindrical shell shall he obtained in the three functions u, v, 
and w. The last of equations (6.2) drops out and the remaining 
ones, in passing to the relative coordinates so that A * B = R, 
may he represented with the aid of table 2. 
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The system of the three equations (8.7) may he reduced to an 
equivalent single equation of the eighth order. Following Galerkin 
(references 7 and 8), the first two equations (8.7) may he substi- 
tuted for X = Y = 0 hy introducing a new function $ = $(a,p) 
and expressing in terms of this function the displacements u, v, 
and w hy the equations 


v 


u 



dadp 2 





2c 


/ b 5 * + g> 5 $ \ 

\Ba^dp da 2 dp^/ 


- (2+U) 


a 5 $ 

da 2 dp 



(8.9) 


= ^1+2 ^ + — 

(3a 4 da 2 dp 2 Sp 4 J 


The last of equations (8.7) assumes the form 
c 2 (v 2 v 2 + 2y 2 + 1)V 2 V% - 2c 2 (l-u) ( - 


da 4 3a 2 Sp 2 , 


V 2 $ +1 


(i.^4.«i2!2| z .o 

da 4 12 Ec 2 


(8.10) 

Equation (8.10) is the fundamental equation of the circular 
shell. In this equation 

^4 


yV _ + 2 

~v 4 + 2 . 2 > 2 
oa da op 


SP' 


( 8 . 11 ) 


For comparison, there are presented the equations obtained on 
the basis of the existing moment theory of Love. These equations, 
given for example in the hook of Timoshenko, may also he repre- 
sented with the aid of table 3. 
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In comparing the equations of table 3 with equations (8.7) 
(given in table 2), it is noted that in equations (8.12) there are 
absent, in the first place, certain terms arising from the moment, 
namely, the terms with c^ = 8^/( 12R^ ) ; in the second place, the 
differential operators (of the second of the third equation and 
third of the second) are asymmetric. The absence of symmetry in 
equations (8.12) is in contradiction to the fundamental theorems 
of the energostatic elastic body. For this reason, the existing 
theory of shells starts from a number of classical problems of the 
mechanics of elastic bodies. 

The previously mentioned defects of equations (8.12) may lead 
to a fundamental error in the problem of the vibration of shells. 
Given any three independent forms of vibration with the correspond- 
ing displacements u, v, and w and applying the method of Gal- 
erkin to equations (8.12), there is obtained for the frequency of 
the vibrations a cubical equation, which being represented in the 
form of a determinant of the third order (corresponding to the 
mechanical significance of the problem) has an asymmetric struc- 
ture. Due to this asymmetry, two of the vibration frequencies for 
arbitrarily chosen forms of u, v, and w may receive Imaginary 
values, a result that is likewise in contradiction to the theory 
of small vibrations of elastic bodies. 

The absence of symmetry in the equations of the moment theory 
of cylindrical shells was noted in previously published papers on 
the theory of shells and thin rods (reference 10) . In these 
papers are given equations of the strength, the stability, and the 
vibrations of shells of composite systems and rods possessing col- 
lateral auxiliary differential operators of the required functions 
of symmetric structure. The recent works on shells (reference 7), 
which improve to a greater or less degree the moment theory, suf- 
fer from the defects pointed out here. 


9. Spherical Shell. Generalization of Equation of Sophle- 
Germain - Lagrange. - In this case. 


k 1 = k 2 = k = 1 /r = constant 

H = k 
L = 0 


(9.1) 
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From equations (5.8), 


'\ 


Aq = 6 + 2kw + v* 

A^ = - V^w - 2k 2 v + 2kv* 


A g = ~ kA ^ 


Xi = x 2 = ° 


(9.2) 




The system (5.4) leads to an equivalent system of the form 

2 ^ 

(X+2^)v 2 0 + 2|ak 2 0 + 2(\-Hi)kV 2 v - (\+2}i) kV 2 V 2 v + \V 2 v* = 

lu 


i -L A (BX) + A (AY) 
5 AB da dg 


2(X-Hx)k0 - (\+2|a) || kV 2 0 + (\+2|a) ^ V 2 Vv 2 + 2n ^ + 

2 

4(\-tp)k 2 v - (3\+2|a) — kV 2 v* + 2\kw* = — Z 

12 8 

2 2 

\0 - (3\+2n) kV 2 w + 2\kw - (i — V 2 v* + (\+2|j.)v* = — Z* 

12 12 8 


V 2 x + k 2 X = 


2|jl8 AB 


t <Br) - k <m) 


j 

(9.3) 


where 9 is the volume dilation of the shell for the tangential 
deformation and 2X is the normal rotation: 


9 = H \t (Bu) + k <A, j 1 

x ■ m [It (Bt) -k (Au il J 


Q 

The symbol V is the differential operator of the second order 
(operator of Beltrami for the sphere): 
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In deriving equations (9.3), 1 + k^5^/l2 ** 1 is assumed 

"because of the smallness of the term k^5^/l2 as compared with 
one. 

The first three equations of (9.3) form a complete system hav- 
ing a symmetric structure with respect to the functions 6, V 2 v, 
and V 2 w*. The fourth equation, independently of the first three, 
determines the normal elongation. 

In the case of a closed spherical shell under the action of 
normal rotation on the inner and outer surfaces for constant (inde- 
pendent of a,P) intensities of these pressures, the differential 
terns drop out, X * Y = 0 and then 

0 - X - 0 

4(A+n)k 2 w + 2\kw* = |- 

2Xkw + (\+2(i)v* « ~ 

o 

where "by virtue of equations (5.5) and (5.7) 

Z = I (1+k?) 2 a 


Z* = | (1+kr) 2 


If w* is set equal to zero and correspondingly the third of 
equations (9.3) is neglected, then for Eu/(1-P 2 ), pi = E/2(l+u) 
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V 2 0 + (l-u)k 2 0 


5^ ,.„2„2 


k7 7 v + (l+u)k7 w = 


1 - 1 ? 1 Tb 


16 AB I da 


(BX) + £ (AY) 


-2 _2 fi 2 l 

- — kV 2 0 + ( 1+ u)k0 + — 7 2 7 2 w + (l-u)— k 2 7 2 w + 2(l-n?)k 2 w = — — 


7^X + k^X = - 


1+v 1 d 


E6 AB da 


(BY) - (AX) 


(9.8) 


For k = 0, the first two of these equations break down into 
the equation of Lame 


7^0 = - 


i-u i a 


18 AB da 


(BY) + (AX) 


(9.9) 


which refers to the problem of the two-dimensional stress state of 
a plate of thickness 8, and the equation of Sophie -Germain- 
Lagrange 

n 2 n 2„ _ -L 17 fr\ - ( 9 . 10 ) 


7 7 w = p Z 


W) 


which refers to the problem of the bending of a plate. 


By eliminating the function 0 from the first two equations 
of (9.8), there is obtained 

— (7 2 7 2 7 2 v + ik^^^w) + k 2 (7 2 w + 2k 2 w) ^ 

12(1-1?) . 


= i J (l+v) ^ hr- (bx) + I- (ay) 


AB da 


(l-Ujl^Z + 7 2 [Z - 


(z - -L A (BX) + A. (AY) 
y 12 AB da dp J /J 


(9.11) 


Equation (9.1l) is the fundamental equation of the spherical 
shell with inextensible normal element and constitutes a natural 
generalization of the equation of Sophie Germain -Lagrange for 
the bending of a plate. 
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Having determined the deflection w and the normal expansion 
from the third of equations (9.8), the tangential displacements u 
and v can he determined h y the equations 

'N 


U = - 


V = - 


1 1 _S_ 

d-„)k 2 ASa 

1 1 _d_ 

(l-u)k 2 B 


( 2 y 

0 + 2kw- — kv 2 w 

v 12 j 


&2 2 ' 
0 + 2kw - — kV w 

i 12 , 


1 SX + 1 1 dw 

^2 B Sp k A da 

_JLldX lldw 

^2 A da k B dp 


l+u 

E8k 2 

1+U 


> 


ESk 


J 

(9.12) 


where 


0 = - 


b c 


12(l+u)k 


V 2 V 2 w - 


S 2 k 
6 ( 1+ u) 


V 2 w 


2kw 


(1-U)5 1 
121 


c 


AB I £ <BX) + I - (AY) 


>1 + i^ 

J E5k 


(9.13) 


The theory of the spherical shell for thick shells (equa- 
tions (9.3)) as well as for moderately thick shells (equations (9.8) 
and (9.11)) have been presented. The fundamental functions chosen 
0, w, X, and w* are invariant relative to the direction of 
the coordinate lines a,p passing through a given point on the 
sphere. 


It follows that the equations given are valid for any system 
of coordinates on the spherical surface. The choice of coordinates 
determines only the differential operator V 2 . If for the coor- 
dinates a, 3 the geographical coordinates were taken, taking a 
as the latitude and 3 as the longitude, then for k = l/R, A = E 
and B = R sin a so that 


V 


2 


E 2 V3a2 + 


cot a ^— + 
oa 


i ^ \ 

sin 2 a dp 2 ' 


(9.14) 


For an arbitrary load (nonsymmetrical problem), equations (9.3) 
and in the particular case (for w* = 0) equations (9.8) or (9.11) 
are integrated by the method of separation of variables with 
respect to the variable p in trigonometric functions and with 
respect to the variable a in Legendre functions . 
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10. General Equations for Stability of Shells - Special Cases. 
It is assumed that the shell has a given system of stresses char- 
acterized only by the tangential (normal and shearing) forces Tp^, 
and S® and in equilibrium with the external forces. It 
shall be considered that the external forces are given with an 
accuracy up to one parameter, for example, the intensity of any 
of the components of the external load. By assigning different 
values to this parameter, different stress states are obtained. 

In a particular case, the internal forces Tp°, T2 and s ° 
may be proportional to the intensity of the external load. For a 
certain value of the load parameter, the equilibrium of the shell 
becomes unstable. 

From the stress state TpO, and S®, the 8 kell passes 

to another state Tp° + Tp, T2 0 + Tg, S® + S, Gp, G2, H, Np, 
and Kg where Tp, T2, . . . , N2 are internal forces arising 
on the loss of stability. It shall be assumed that the forces Tp, 
T2, . . . , N 2 and the corresponding deformations are infinitely 
small magnitudes. Because the change in the deformed state of the 
shell, associated with loss in stability, is characterized by a 
change in form of the middle surface, it is necessary, in order to 
obtain the equations of stability, to take into account the varia- 
tions of the magnitudes referring only to the second -quadratic form 
of the surface. 

The stability equation is obtained from the equations (6.2) 
and (6.3) given for shells of medium thickness or from equa- 
tions (7.4) for thin shells. It is necessary in the first place 
to refer all static and kinematic magnitudes entering these equa- 
tions to variations of the stress and the deformed state of the 
shell that arise on loss of stability and in the second place to 
consider the components X, Y, and Z as those surface forces 
that are obtained when an element of the shell AB da dp with the 
contour forces Tp° + Tp, T2 0 + T2, S° + S, . . . is carried 
into the new deformed state determined by the displacements u, v, 
and w. 


With the passage of the shell into the deformed state, the 
normal to the middle surface will have a new direction determined 
by the angles of rotation 


<11 


(*2 v “ I |r) 




J 


q 2 = - 


(10.1) 
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relative to the tangents to the lines a » constant and p = constant 
of the initial state. For the components X, Y, and Z of the 
vector of the reduced surface force on the axes of coordinates a, 
p , and 7 of the movable trihedron of the middle surface, small 
magnitudes are readily obtained with an accuracy up to second order 



By substituting the values of X, Y, and Z thus obtained 
in equations (6.2) and neglecting in these equations and m^, 

the general equation of stability of the shells is obtained. In 
the case of a thin shell, X, Y, and Z must be substituted in 
equations (7.4). 

In either case, there is obtained with the accuracy of the 
load parameter a complete system of homogeneous differential equa- 
tions in the required functions that determine the deformations of 
the shell associated with loss in stability. To this system are 
added the boundary conditions (homogeneous). 

The critical stress state Tj°, T 2 °, and S° determined by 
the parameter of the external load entering linearly in equa- 
tions (10.2) is thus determined by solving the homogeneous boundary 
problem described here. 

Inasmuch as of the three variables u, v, and w the normal 
displacement w has the principal effect on the change in shape of 
the shell, in equations (10.2) the terms with the tangential dis- 
placements u and v may be neglected. Then 
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(10.3) 


In the case of a thin shell, the components X and Y, being 
proportional to the curvatures k-j_ and kg, may be taken equal to 
zero. On the basis of equations (7.8) 


~ V e 2 V e 2 $ - (HV e 2 -LV h 2 )v = 0 



(10.4) 

These equations constitute the general equations of stability 
of thin shells in the two functions $ and w and permit deter- 
mining the critical stresses for very general assumptions with 
regard to the given stress state. 

The general theory of the stability of shells has been pre- 
sented exactly as given by equations (10.2), (6.2), and (6.3) and 
approximately as given by equations (10.3), (7.4) or (10.4). 

This theory represents a considerable generalization of a num- 
ber of problems on stability of elastic systems, starting with the 
simplest problem of longitudinal bending and ending with the sta- 
bility of shells of arbitrary shape for arbitrarily given initial 

stress state Tq° = T]°(a,p), T 2 ° = T 2 °(a,g), and S° * S°(a,p), 
the critical value of which is determined. 
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Thus, for example, in equations (10.4) setting k^ = kg = 0, 

from the second equation (the operator 

2 2 ki+ko 2 k-i-ko 2 . 

HV^ - LV>, = — - - V- - _± — £ V, becomes zero), the 

e - n 2 ® 2 

equation of stability of a plate in arbitrary orthogonal curvilin- 
ear coordinates a and 8 are obtained. For 



there is obtained 

^ - T ,° T,° . 2S° » 0 (10.5) 

12(l-U 2 ) 5a 58 5a58 


Equation (10.5) is the well-known equation, in rectilinear 
coordinates, of the stability of a plate loaded by forces on the 
boundary. 

For h-j_ = kg = k = l/R = constant, S® = 0, and 
T-j® = Tg® = -pR/2 = constant, the following equation is obtained: 


R 2 5 2 

12(l-u 2 ) 


2 2 2 
7 e ^e 7 e w + 


2 

V e w + 


pR 3 

2E5 


2 

7 e w 


0 


( 10 . 6 ) 


which refers to the stability of a spherical shell of radius R, 
under an internal pressure p = constant, the parameter of this 
equation being p. 


A cylindrical shell will be considered, starting at first 
from the more accurate equations (10.2) and (8.7). If a and 6 
are the absolute coordinates, A = B «= 1. Equations (10.2) for 
k^ a 0 and kg = constant become 
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(10.7) 

ST] 0 iqO ST ? ° £q0 

If — — + £2- = 0, -r-=- + ~ =0, and A = B » E (hence 

5a Bp ' dp da ' 

a and p are relative coordinates ) , equations (10.7) become the 
following: 



On the basis of equations (10.8), table 2 (equations (8.7)) 
assumes the form given by table 4. Equations (10.9) (table 4) in 
the secondary operators possess also in this case a symmetrical 
structure, a fact that as already noted is in agreement with the 
theorem of reciprocity and therefore the critical forces will 
always be real.-*- 


%he equations that are used by Timoshenko (reference 3) and 
other authors (references 7 and 8) are assymetric with respect to 
the secondary terms and consequently do not correspond to the funda- 
mental theorems of the theory of elasticity. 
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These equations are the general equations of stability of a 
circular shell, obtained in strict correspondence with the funda- 
mental hypothesis of Kirchhoff - Love and make it possible to con- 
sider a number of problems of practical interest on the determina- 
tion of the critical loads of the shell. 


Because the components X, Y, and Z determined in the gen- 
eral case of the stress state by equations (10.8) are obtained 
with account taken of the exact values of the angles of rota- 
tion q^ and q 2 of the normal of the shell, equations (10.9) 
are applicable also to shells of medium thickness . 


In the case of a shell for which 5/R $ l/ 3®, the tangential 
contour force Y represents, according to equations (10.7), a mag- 
nitude that is small in comparison with the normal force Z. The 
tangential contour displacement v on deformation of the shell 
accompanied by the change in shape of the cross section is a mag- 
nitude that is likewise small compared with the normal displace- 
ment w. By assuming for a thin shell the magnitude Y to be 
equal to zero and neglecting in the last of equations (10.8) the 
tangential displacement v, 


z = T 

r z 


“_d 

f m 0 dw 

o 0 dw"\ 

d 

f m 0 dw 

„0 dw\ 

da 


+ s 

+ dP 


+ 3 sJJ 


( 10 . 10 ) 


The general stability equation of a shell for given assumptions 
as to the force Y may be obtained from equation (8.10) by sub- 
stituting in this equation the value of Z determined by equa- 
tion (10. 10). This equation has the form Z 

c 2 (vW + 1 )vH - 2cHl-v)( *2# + (i-p2) ** - 

Vda da dp / da L 


- JL (t-i 0 — v^+ s° — + — f Tp° — v 4 ® + s° A. = 0 

EB Ida \ da dp / dp \ dp da /J 


( 10 . 11 ) 


p 

If in equations (10.11) the second and third terms of the 
first component and the complete second component are neglected, 
the approximate equation for the stability of a thin cylindrical 
shell shall be obtained. 
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Equation (10.11) for these assumptions is equivalent to the 
system of equations (10.9) in the three functions. The displace- 
ments u, v, and w are determined in terms of the fundamental 
function $ by equation (8.9). It should be noted that 
T 2 ° = S° = 0, that is, in the case of the stress state character- 
ized only by longitudinal normal forces T]® (central compression, 
for example, pure bending, eccentric action of longitudinal com- 
pressive or tensile forces, and so forth), these assumptions drop 
out. The equations given here are the general equations of the 
stability of a cylindrical shell from which the critical stress 
can be determined for very different assumptions both as regards 
the given external forces and as regards the boundary conditions. 
Thus, for example, the equations of stability can be obtained for 
the f olio-wing cases : 

1. Central compression of a shell by a force P 



2. Pure torsion of a closed shell 



3. Shells under the action of an external normal pressure and 
immovably clamped at the longitudinal edge against displacements u 
and w 



4. Shells under the simultaneous action of a longitudinal com- 
pressive (or tensile) force P and twisting moment M 
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In this case for the load parameter, there may "be taken the 
magnitude P for a given value of the magnitude M or M for a 
given value of P, or the ratio of these values as a function of 
the conditions of the problem 

5. Shells under the action of only a single bending moment 
(pure bending) or of a moment and a longitudinal force (bending 
with tension or compression) 

6. Shells under the action of a transverse load producing at 
the sections a = constant longitudinal, normal, and shearing 
forces T-j 0 and S°, determined by the usual elementary theory 

of the bending of beams, and so forth 

In all of these cases except cases 5 and 6 the differential 
equations of stability have constant coefficients . 

The critical stresses are determined by solving the homogen- 
eous boundary problem by equation (10.11) or in the case of a more 
accurate solution by the system of equations (10.9) and the homo- 
geneous boundary conditions . If the shell of length Z on each 
of the curvilinear edges a « 0 and a = Z/R is hinge -supported 
on a diaphragm that is rigid in its plane and flexible in the trans- 
verse plane, the function $ corresponding to these boundary con- 
ditions may be approximated for the closed shell in the form of a 
double trigonometric series : 

$ = A mn sin cos n3 (m,n = 1, 2, 3, . . .) 

and for shells of open profile in the form of a trigonometric series 
in only one variable a: 

® = V * m (P) sin 52®SL (m . 1, 2, 3, . . .) 

^ x Z 


where the function ^(0 ) is determined by ordinary differential 
equations (homogeneous with one parameter) and the boundary con- 
dition (likewise homogeneous), which must be given on the straight 
edges of the shell- 
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REMARKS 

The theory given ia of general character and permits solving 
a number of practically important problems on the strength, the 
stability, and the vibrations of shells. Thus for example: 

1. Computation of shells by the method of the theory of com- 
plex variables . - In reference 12, it is shown that for shells char- 
acterized by middle surfaces of the second order vith positive 
Gaussian curvature (spherical, elliptical, and parabolical) the 
equations of the momentless theory characterized by the mixed dif- 
ferential operator HV e 2 - Lv^ 2 leads, through transformation of 
the independent variables, to the Cauchy -Re iemann equations. 

These investigations show that the more accurate equations (7.8) 
relative to the moment theory of thin shells will be of the ellip- 
tic type for middle surfaces of the second order. These equations 
for such surfaces also lead to the equations of Cauchy -R iemann. 

It then follows that the computation of such shells by the moment 
theory may be effected by the methods of the theory of functions 
of a complex variable by developing and generalizing the known 
methods of Muskhelishvili (reference 13) on the two-dimensional 
problem of the theory of elasticity. In particular, it is of 
interest to determine the stresses and the deformations of shells 
of spherical, elliptic, and parabolic types due to the action of 
a concentrated force applied at any point of the middle surface. 

2, Circular cylindrical shell under the action of a concen- 
trated force . - The solution of this problem may be obtained by the 
integration of equation (8.10) or for the thin shell of equa- 
tions (7.8) by the method of separation of variables. (In this case 

HV e 2 - L7 h 2 = *2 $—.) The functions required me, be approximated 

by trigonometric series in one of the variables a or p, as in 
the method of Failon-Ribier for the two-dimensional stress state of 
a rectangular plate and in the method of Morris -Levy for the case 
of the bending of such a plate. The Green Function may be repre- 
sented by a Fourier Integral. 

5. Tension in a closed circular shell having somewhere on the 
surface an opening of given shape . 

4. Torsion of a circular shell weakened by an opening . - Both 
of these problems may be solved also with the aid of trigonometric 
series . 
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5. Stability of an open circular shell in the oaae; (a) cen- 
tral compression, (b) pure bending, (o) compression with bending, 
and (d) bending by given transverse forces . - This problem la 
solved by applying to the stability equations given ordinary trig- 
onometric series in the variable along the generator. 

6. Stability of a closed circular shell in torsion . - The 
required functions in this case can be given in the form of trig- 
onometric series in the variable p (in the direction of the trans- 
verse arc). 

7. Stability of a spherical shell under the action of an external 
hydrostatic pressure . - The differential equation corresponding to 
this problem can be integrated by the method of separation of the 
variables by applying trigonometric functions and functions of 
Legendre . 


Translated by S. Eeiss 
National Advisory Committee 
for Aeronautics. 
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